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Abstract 

The theory of plasma physics offers a number of nontrivial examples 
of partial differential equations, which can be successfully treated with 
symmetry methods. We propose the Grad-Shafranov equation which 
may illustrate the reciprocal advantage of this interaction between plasma 
physics and symmetry techniques. A symmetry classification of the Grad- 
Shafranov equation with two arbitrary functions F (u) and G(u) of the 
unknown variable u = u(x, t) is given. The optimal system of one- 
dimensional subalgebras is performed. This latter provides a process for 
building new solutions for the equation. 
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1 Introdutction 

The equation we are going to examine contains indeed two arbitrary functions 
F(u), G(u) of unknown variable u = u(x, t), and the goal is now to perform the 
symmetry classification of this equation, i.e. to find those F, G for which the 
equation admits nontrivial symmetries. In general, the symmetry propoerties 
of an equation may strongly depend on the choice of the arbitrary functions 
involved. 

The PDE we want to consider is 

u xx + a-u x + u tt = x p F(u) + G(u) (1) 

X 

with a = — 1 , b = 1 this equation is known in plasma physics as the Grad- 
Shafranov equation^ see [?]) and describes the magnetohydrodynamic force bal- 
ance in a confined toroidal plasma. In this context, u is the so-called magnetic 
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flux variable, a; is a radial variable, then x > 0, while the two arbitrary flux 
functions F(u),G(u) are related to the plasma pressure and current density 
profiles. 

This equation admits well known properties and a rich literature is devoted 
to it. After some comments on the general peculiarities of the equation and 
of its symmetry properties, we shall provide the algebra of its exact Lie point 
symmetries and probe some properties of its algebra. Next, we shall study 
one-dimensional optimal system of subalgebras. 

According to the standard definitions and procedure (see [?]), we can now look 
for the equivalence group; preliminary, we look for the full groups of the equation 
([2]), i.e. the intersection of all groups admitted by ([2]) for any arbitrary choice 
of F and G. It turns out that this kernel is almost trivial: it contains indeed 
only the transformation of the variable t. 

First of all, we exclude from our classification the case a — and p — (or 
a = and F(u) = 0) because in this case our equation coincides with well 
known nonlinear Laplace equation. 

Now, we investigate the case, F = 1 and G = for the Grad-Shafranov equation. 



2 Symmetry methods 

Let a partial differential equation consists p independent and q dependent vari- 
ables. The one-parameter Lie group of transformations 

x t = Xi + sC(x, u) + 0(s 2 ); u a — u a + SLp a (x,u) + 0(s 2 ), (2) 

where £' = ^\ s=0 ,i = l,---,p, and ip a = ^±f\ s=0 ,a = 1, • • • ,q, are given. The 
action of the Lie group can be recovered from that of its associated infinitesimal 
generators. We consider general vector field 

i— 1 a— 1 

on the space of independent and dependent variables. Therefore, the character- 
istic of the vector field Y given by © is the function 

p du a 



Q a (xM 1) ) = V a (x,u)-J2eM-^> ol = 1,---, q. (4) 
The second prolongation of the infinitesimal generator 



d d d 

X = £\x,t,u)— + £ 2 (x,t,u)—+ ip(x,t,u)— (5) 

is the following vector field 



du x du t du xx du xt du u 
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with coefficients 





= D x Q-\ 


- £ u xx + i 2 u xtl 
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= D t Q 4 


■ ^Uxt + £ 2 u tt , 


lf xx 


= D 2 X Q A 


£ Uxxx 4 £ ^xxt: 


tp xt 


= D x D t C 


3 + £}u xxt + Sj, 2 Uxtt 


v u 


= A 2 Q^ 


-£ l u xtt + £, 2 u tt t, 



(6) 



where the operators D x and £>t denote the total derivatives with respect to x 
and t : 

By theorem 6.5. in [?], the vector field X is a one parameter of Grad-Shafranov 
equation if and only if 

(2) 1 1 

Pr X(u xx u x + Uu — x 2 ) — whenever u xx u x + u u — x 2 — 0. (7) 

x x 

So, we apply the criterion of infinitesimal invariance in order to determine sym- 
metries of the Grad-Shafranov equation[?]. Therefore, the infinitesimal symme- 
try criterion is 

^ - ^u x - V + ip u - 2x£ 1 - (8) 
ar x 

Substituting the formulas (J?) into ©, we are left with a polynomial equation 
involving the various derivatives of u whose coefficients are certain derivatives 
of £ , £ 2 and <p only depend on x,t,u we can equate the individual coefficients 
to zero, leading to the complete set of defining equations: 

2 2 H 2 2 & 

£u = Oj^tt — ? Ctx = 0? Cxa; = i (9) 

a; X 

4£f x 3 + <p x - ip u x 3 - (p xx x 
<Ptt = , (10) 

X 

ptu = —zr, <fuu = Q, ¥ux = 0, £ (x,t,u) =£ t x. (11) 

Lx 

By solving this system of PDEs, we state the following theorem. 

Theorem The Lie algebra q of the symmetry group G associated to the Grad- 
Shafranov equation is generated by the vector fields 

Xi = d t , 

x 4 

X 2 = xd x + td t + -^dm 

X 3 = (~+u)d a , (12) 

X 4 = txd x + ht 2 -x 2 )d t + ^-t(7x 4 + 8u)d u , 
I lb 

X 5 = ip(x,t)d u . 
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The commutation relations of the 4-dimensional Lie algebra g spanned by the 
vector fields Xi,X 2 , X 3 , X4 are shown in the following table. 



Table 1: Commutation relations satisfied by infinitesimal generators in (33). 
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x 2 


x 3 


x 4 


x 1 





x 1 





X 2 + ^X 3 


x 2 


-x 1 








Xi 


x 3 














Xi 


—X 2 — ^x 3 


-x A 









3 The Lie algebra of symmetries 

In this section, we determine the structure of full symmetry Lie algebra g of 
Eq. ©. 

Theorem The full symmetry Lie algebra g of Eq.^j has the following semidi- 
rect decomposition: 

g = R x gi (13) 
where gi is a semi-simple Lie algebra. 

Proof: The center z of g is Spann{X 3 }. Therefore the quotient algebra gi = 
g/z is Span^{Yi,Y 2 ,Y 3 }, where Yi — Xi + z for i = 1,2,3. The commutator 
table of this quotient algebra is given in the following table: 



Table 2: Commutation relations satisfied by infinitesimal generators in (33). 



[,] 


Y t 


Y 2 


Y 3 


Yi 





Yi 


Y 2 


Y 2 


-Yi 





Y 3 


Y 3 


-Y 2 


-Y 3 






The Lie algebra g is non-solvable, because 

g (1) = [g,g] =Span n {Xi,X 2 + ±X 3) X 4 }, (14) 

g (2) = [g (1) ,g (1) ]-g (1) . ' (15) 

Similarly, gi is semi-simple and non-solvable, because 

s{ 1} = fei.gi] = Spann{Yi,Y 2 ,Y 3 } = gl . (16) 

The Lie algebra g admits a Levi decomposition as the following semi-direct 
product g = r K s, where r = Spann{X 3 } is the radical of g (the largest 
solvable ideal contained in g), and s = Spann{Xi,X 2 + ^X 3 ,Xi}. 

The ideal r is a one-dimensional subalgebra of g, therefore it is isomorphic to 
R; Thus the identity g = r K s reduces to g = R x gi . 
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4 Optimal system of subalgebras 



Consider a system of partial differential equations A defined over an open subset 
M C X x U ~ MP x M. q of the space of independent and dependent variables. 
Let G be a local group of transformations acting on M. Roughly, a solution 
u = f(x) of the system is said to be G- invariant if it is left unchanged by all the 
group transformations in G. If G is a symmetry group of a system of partial 
differential equations A, then, under some additional regularity assumptions on 
the action of G, we can find all the G-invariant solutions to A by solving a 
reduced system of differential equations, denoted by A/G, which will involve 
fewer independent variables than the original system A. In general, to each s- 
parameter subgroup H of the full symmetry group G of a system of differential 
equations in p > s independent variables, there will correspond a family of 
group-invariant solutions to the system. Since there are almost always an infinite 
number of such subgroups, it is not usually feasible to list all possible group- 
invariant solutions to the system. We need an effective, systematic means of 
classifying these solutions, leading to an "optimal system" of group-invariant 
solutions from which every other such solution can be derived. Since elements 
5 £ G not in the subgroup H will transform an iJ-invariant solution to some 
other group-invariant solution, only those solution, not so related need be listed 
in our optimal system. 

Let G be a Lie group. An optimal system of s-parameter subgroups is a list 
of conjugacy inequivalent s-parameter subgroups with the property that any 
other subgroup is conjugate to precisely one subgroup in the list. Similarly, a 
list of s-parameter subalgebras forms an optimal system if every s-parameter 
subalgebra of is equivalent to a unique member of the list under some element 
of the adjoint representation: I) = Adg(t)),g € G. 

Proposition 3.7 of [?] says that the problem of finding an optimal system of 
subgroups is equivalent to that of finding an optimal system of subalgebras. 
For one-dimensional subalgebras, this classification problem is essentially the 
same as the problem of classifying the orbits of the adjoint representation, since 
each one-dimensional subalgebra is determined by a nonzero vector in g. This 
problem is attacked by the naive approach of taking a general element V in 
q and subjecting it to various adjoint transformations so as to "simplify" it as 
much as possible. Thus we will deal with the construction of the optimal system 
of subalgebras of q. 

The adjoint action is given by the Lie series 

Ad(exp(eF i )F J ) = Y j - e[Y u Yj] + e 2 [Y u [Y u Yj]] - • • • , (17) 

where [Yj, Yj] is the commutator for the Lie algebra, e is a parameter and i, j = 
1,2,3. 

The adjoint representation of q is listed in the following table, it consists the 
separate adjoint actions of each element of g on all other elements, 
where the (i,j)-th entry indicating Ad(exp(eYi)Yj) . 
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Table 3: Adjoint relations satisfied by infinitesimal generators in (33). 




A4 X\ — eX2^sX^ H~ ■^6^X4 X2 — eX\ A3 X4 



Theorem An optimal system of one dimensional Lie subalgebras of Grad- 
Shafranov equation is provided by those generated by 



where a, b G R are arbitrary constants. 

Proof Let g is the symmetry group of Eq. ([3]) with adjoint representation 
determined in Table 2 and 



is a nonzero vector field of g. We will simplify as many of the coefficients of 
a,, i = 1 , • • • 4; as possible through judicious applications of adjoint maps to X. 



Suppose first that 04 ^ 0. Scaling X if necessary, we can assume that 04 = 1. 
Referring to table 2, if we act on such a X by Ad(exp(— (I3A1)), we can make 
the coefficient of A3 vanish. Thus, every one-dimensional subalgebra generated 
by a X with 04 ^ is equivalent to the subalgebra spanned by aX\ + 6A2 + X4 
where a, b £ R are arbitrary constants. No further simplifications are possible. 



The remaining one-dimensional subalgebras are spanned by vectors of the above 
form with 04 = 0. If 03 =/= 0, we can scale to make 03 = I. There are two 
subcases. 

Case 2.1: 

If 02 7^ 0, then we can cancel the coefficient of X\ by acting on X by Ad(exp(— ^X\)). 
So that, X is equivalent to a scalar multiple 0X2 + A3 for some a € R. 

Case 2.2: 

If a2 = 0, we can act by adjoint map generated by X2 to arrange the coefficient 
of Xi either +1, -1 or 0. Therefore, any one-dimensional subalgebra spanned 
by X with 04 = 0, 03 = 1, 02 — is equivalent to one spanned by either A3 + 
Ai, A 3 — Ai or A 3 . 



1) X 1 , 2)A 2 , 3)A 3 , 4)A 3 -A 1; 5)A 3 + X 1 
6) aA 2 + A 3 , 7)aX 1 + bX 2 + A 4 . 



(18) 
(19) 



A = aiXi + 02X2 + &3A3 + a 4 A 4 



(20) 



Case 1: 



Case 2: 



Case 3: 
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The remaining cases, = = 0, are similarly seen to be equivalent cither to 
-^2(02 7^ 0) or to X\{a 2 = a 3 = 04 = 0). 

There is not any more possible case for studying and the proof is complete. 

In continuation, we find some group invariant solutions of the equation(??) 
corresponding to 1-dimensional subalgebras generated by Xi,X 2} X 4 . 



Consider 



invariants are 



x 2 = xd x + td t + — d u (21) 



Ci = -, C 2 = ~lx 4 + u (22) 
x 8 



So, group invariant solution associated to these invariants has the form 



(t + Cl)vFT?+Cat m 



In the case of X4, invariants are 



t 2 +x 2 8u-x 4 



therefore, we get the associated group invariant solution 



\{xH 2 + x 7 )y^±Z+ (2C 2 t 2 + d)x 5 / 2 + C 2 (x^ 2 + y^t 4 ) 

^±Zx(t 2 +x2) 

In the case of X = X 2 + X 3 , we obtain the following invariants 



Thus corresponding invariant solution is 



x 4 „ It 2 + x 2 



(25) 



u {x,t) = — + c 2 y —^—x + at (27) 

In the case of X = X\ + X3, we get the invariants 

a = x, C 2 = — -(x 4 - 8it) exp(-i) (28) 
8 

therefore, we obtain the following invariant solution associated to these invari- 
ants. 

u(x, t) = x(^-x 3 + a exp(i)BesselJ(l, x) + C 2 exp(t)BesselY(l, a;)) (29) 



where BesselJ and BesselY are the Bessel functions of the first and second kinds, 
respectively. They satisfy Bessel's equation. 

For the vector field X\ — d t , global invariants are C\ = x, C2 — u. So, the 
solution of reduced equation is 

u = x 4 + AC x x 2 + C 2 (30) 

5 New admitted symmetries 

In this section, we probe the remainder cases. 

If F = exp(2u) and G — exp(u), there is a new admitted symmetry which has 
the form 

X = xd x + td t - 2d u (31) 

If F = and G = u 1+ « for all q ^ 0, there is a new admitted symmetry 

which is 

X = xd x + td t - 2qud u (32) 

If F = 1, G = u, new symmetries are 

X = (x 2 + u)d u , X = if>(x, t)d u (33) 

Furthermore, the associated invariant solution for the global invariants C\ = 
x, C 2 = u is 

u(x, t) = -x(x - C 2 BcssclI(l, x) + CiBessclK(l, x)) (34) 

where Bessell and BesselK are the modified Bessel functions of the first and 
second kinds, respectively which satisfy the modified Bessel equation. 

If F = u, G = 1, the new symmetry has the form 

X = (u + i/>(x,t))d u (35) 

In addition, the corresponding invariant solution for the global invariants C\ = 
x, C2 = u is 

u(x,t) = cosh(^)(C 1 - ^Shi(^)) + sinh(^)(C 2 + ^Chi(^)) (36) 

where Shi and Chi are hyperbolic Sine integral and hyperbolic cosine integral, 
respectively. 



Shi(x) 



r^dt (37) 

Jo t 

Chi(x) = 7 + ln(.x)+ / X cosh (*)-l dt (38) 
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